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SUMMARY 


The  frequency  behavior  of  the  scattered  acoustic  field  produced  by  a  plane  wave 
impinging  on  an  elastic  body  immersed  in  an  infinite  fluid  medium  is  dominated  in 
certain  frequency  ranges  by  large  peaks  arising  from  resonant  modes  in  the  elastic 
structure.  Computational  models  have  been  developed  at  the  Naval  Ocean  Systems 
Center  for  solving  these  elastic  scattering  problems.  This  report  shows  that  the 
scattered  pressure  can  be  represented  as  a  sum  of  contributions  from  in  vacuo  modes 
of  the  structure.  The  spectral  response  of  individual  terms  in  this  summation  can  be 
used  to  identify  which  modes  contribute  to  each  peak  in  the  overall  spectral  response 
of  the  scattered  field.  In  many  cases  a  small  number  of  modes  dominate  the  response 
in  the  vicinity  of  these  peaks. 
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INTRODUCTION 


The  frequency  behavior  of  the  scattered  acoustic  field  produced  by  a  plane  wave 
impinging  on  an  elastic  body  immersed  in  an  infinite  fluid  medium  is  dominated 
in  certain  frequency  ranges  by  large  peaks  arising  from  resonant  modes  in  the  elas¬ 
tic  structure.  A  technique  for  computing  the  scattered  acoustic  field  was  described 
in  a  previous  Naval  Ocean  Systems  Center  (NOSC)  technical  report  (Schenck  and 
Benthien,  1989).  This  technique  makes  use  of  a  finite  element  model  of  the  elastic 
structure  and  a  Helmholtz  integral  model  (CHIEF)  (Benthien,  Barach,  and  Gillette, 
1988)  (Schenck,  1968)  of  the  acoustic  medium.  The  purpose  of  the  present  report 
is  to  show  how  quantities  computed  with  the  above  technique  can  be  used  to  iden¬ 
tify  which  in  vacuo  structural  modes  are  the  main  contributors  to  each  peak  in  the 
scattered  field  frequency  response. 

The  second  section  of  tliis  report  contains  a  brief  survey  of  the  numerical  approach 
used  to  solve  the  acoustic  scattering  problem  as  well  as  the  basic  equations  needed 
for  identifying  modal  contributions.  The  third  section  contains  numerical  results 
obtained  by  applying  the  modal  identification  technique  to  two  scattering  problems. 
The  fourth  section  contains  a  brief  summary  of  conclusions  which  illustrate  the  utility 
of  this  technique  in  understanding  the  underlying  structure  of  the  scattered  field. 
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BASIC  EQUATIONS 


Throughout  this  section  all  pressures,  displacements,  and  velocities  will  be  repre¬ 
sented  by  iheir  complex,  frequency  dependent  Fourier  components  corresponding  to 
the  time  dependence  The  elastic  structure  is  modeled  using  the  finite  element 
equations  (Zienkiewicz,  1977) 

i-ioHl  +  K)U  =  F,  (1) 

where  M  is  the  mass  matrix,  K  is  the  stiffness  matrix,  F  is  the  load  vector,  and 
is  a  vector  whose  components  are  the  displacement  degrees-of-freedom  at  all  the 
nodes  in  the  body.  In  the  scattering  problem  under  investigation,  the  load  vector  F 
is  completely  determined  by  the  acoustic  pressure  p  on  the  wet  surface  S  of  the  body. 
The  components  of  F  are  given  by 

Fm  =  -  f^p(t>^-ndS,  (2) 

where  p  is  the  acoustic  pressure,  n  is  a  unit  normal  to  S  pointing  into  the  fluid,  and 
4>m  is  a  finite  element  vector  interpolation  function.  The  displacement  u{x)  at  a  point 
X  in  the  body  can  be  expressed  in  terms  of  the  interpolation  functions  as  follows: 

M 

=  Y.  (3) 

m=l 

where  Ui,. . .  ,Um  are  the  components  of  U . 

The  CHIEF  formulation  of  the  acoustic  part  of  the  problem  is  based  on  the 
Helmholtz  integral  relations.  The  surface  Helmholtz  integral  equation  is  approxi¬ 
mated  by  the  system  of  algebraic  equations 

AP  =  BV  ^  (4) 

where  the  matrices  A  and  B  involve  integrals  of  the  free-space  Green’s  function  and 
its  normal  derivative  over  the  subdivisions  5„  of  5;  P  and  V  are  vectors  whose  nth 
components  are  the  pressure  and  normal  velocity  (assumed  to  be  constant)  on  the 
subdivision  5„  of  5';  and  E'"'’  is  a  vector  (Benthien,  et  ah,  1988)  whose  nth  component 
is  the  value  of  the  incident  pressure  wave  at  a  reference  point  on  S^-  The  scattered 
pressure  p’  at  a  field  point  x  exterior  to  the  body  can  be  approximated  by 

p%x)  —  a^{x)P b^{x)V^  (5) 

where  a{x)  and  b(x}  are  vectors  whose  components  (Benthien,  et  ah,  1988)  involve 
integrals  of  the  free-space  Green’s  function  and  its  normal  derivative  over  the  subdi¬ 
visions  Sn  and  evaluated  at  the  field  point  x.  Equations  (4)  and  (5)  can  be  combined 
to  give 

;/(x)  =  /(a:)V-hp”(x),  (6) 
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where 


[x]  ~  a^{x)A  ^  B  +  6^(i) 


p”(x)  =  J{x)A-^P'^^. 

It  can  be  seen  from  equation  (6)  that  p"(a;)  is  the  rigid  scattering  from  the  body  (i.e., 
the  scattered  pressure  when  v  is  constrained  to  be  zero). 

Applying  the  assumption  that  p  is  piecewise  constant  on  the  subdivisions  of  S  to 
equation  (2),  gives 

F  =  ~CDP,  (7) 


where 


Cmn  —  Q  f  ^ 

Ofi  JSn 


D  -  diag(5'i,52,...,5N). 

Equations  (1)  and  (7)  can  be  combined  to  give 


U  =  +  K)-^CDP 


Since  different  interpolation  schemes  are  used  in  the  finite  element  model  of  the 
structure  and  the  CHIEF  model  of  the  acoustic  medium,  it  is  impossible  to  enforce 
exact  continuity  of  normal  displacement  across  S.  However,  this  continuity  is  ap¬ 
proximately  enforced  by  equating  the  CHIEF  normal  velocity  Vn  to  the  average  of 
the  finite-element  normal  velocity  over  5„,  i.e.. 


On  JS„ 


iijju  ■  n  dS. 


Combination  of  equations  (3),  (8),  and  (9)  gives 

V  =  luC'^U 

=  +  K)-^CDP. 


Let  E  be  a  matrix  whose  columns  are  the  in  vacuo  normal  modes  of  the  structure. 


I<E  =  MEn,  (11) 

where  =  diag(u;j , . . .  ,lj^)  is  the  diagonal  matrix  of  eigenfrequencies.  Since  the 
modes  are  M-orthogonal,  they  can  be  normalized  so  that 

E'^ME  =  I.  fl2l 
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It  is  easily  verified  that  the  inverse  matrix  -f  K)  ^  can  be  expressed  in  terms 

of  the  normal  mode  matrix  E  as  follows: 

(13) 

This  form  is  convenient  since  the  diagonal  matrix  (— u;^/  +  Q)  is  easily  inverted  at  all 
frequencies.  Combination  of  equations  (10)  and  (13)  gives 

V  =  -iuC^Ei-cj^I  +  ny'^E'^CDP.  (14) 

Equations  (4)  and  (14)  can  now  be  combined  to  give 

[A  +  iuBC'^E{-<Jl  +  fl)-*A^CA)]P  =  (15) 

Equation  (15)  can  be  solved  for  the  surface  pressure  vector  P.  Once  P  is  determined, 
the  velocity  vector  V  can  be  obtained  from  equation  (14)  and  the  scattered  pressure 
p^(x)  can  be  computed  at  any  exterior  field  point  x  using  equations  (5)  or  (6). 


The  remainder  of  this  section  will  be  devoted  to  the  development  of  a  modal 
expansion  for  p'’(x).  Each  column  of  E  represents  an  in  vacuo  mode  of  the  structure. 
If  Cn  is  the  nth  column  of  E,  then  equation  (13)  can  be  written  in  the  alternate  form 

M  1 

+  K)-'  =  ■£  (16) 

m=l  '^m  ^ 

Combination  of  equations  (8)  and  (16)  gives 


In  view  of  equation  (10),  it  follows  that 


Substitution  of  equation  (18)  into  equation  (6)  gives 


p’{x)  -  p”(x 


(19) 


The  nth  term  of  the  sum  in  equation  (19)  is  the  contribution  of  the  mode  e„  to 
the  scattered  field.  Once  equation  (15)  is  solved  for  P,  the  contributions  of  various 
modes  can  be  determined  by  plotting  individual  terms  of  equation  (19)  or  partial 
sums  of  equation  (19)  versus  frequency  and  comparing  the  results  with  the  plots  of 
p^{x)  —  p’'*(x)  versus  frequency. 


4 


NUMERICAL  RESULTS 


In  this  section,  the  results  of  the  previous  section  will  be  applied  to  two  scatter¬ 
ing  problems.  The  first  problem  consists  of  determining  the  scattered  acoustic  field 
produced  when  a  plane  incident  wave  strikes  a  hollow  spherical  shell.  The  ratio  of 
the  shell  thickness  to  the  mean  radius  is  0.03  .  The  material  parameters  of  the  shell 
are 


•  Density  =  7669  Kg/m^ 

•  Young’s  modulus  =  2.07  x  10" 

•  Poisson’s  ratio  =  0.3  . 

The  material  properties  of  the  water  are 

•  Density  =  998  Kg/m^ 

•  Sound  Speed  =  1486  m/s. 


Figure  1  shows  the  contributions  of  the  second  (the  so-called  accordion)  mode  to  the 
magnitude  of  the  form  function.  The  dotted  curve  is  the  sum  of  the  rigid  scattering 
form  function  and  the  form  function  obtained  by  using  the  m  =  2  term  of  equation 
(19).  It  should  be  noted  that  the  eigenfrequency  of  mode  2  in  vacuo  occurs  at  a  ka 
of  2.55,  wherecis  its  contribution  to  the  overall  response  occurs  at  a  ka  of  about  1.58 
due  to  the  mass-like  loading  of  the  acoustic  field.  Figure  2  shows  a  similar  result  for 
the  fourth  structural  mode.  In  fact,  each  of  the  peaks  shown  in  the  overall  response 
is  primarily  due  to  a  single  (higher  frequency)  in  vacuo  mode  of  the  spherical  shell. 
This  simple  situation  is  not  usually  observed  in  the  scattered  field  produced  by  other 
shapes  where  peaks  are  often  due  to  the  interaction  of  several  structural  modes.  This 
will  be  illustrated  in  the  second  example. 
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BACKSCATTERING  FROM  SPHERICAL  SHELL 


Figure  1.  Contribution  of  mode  2  to  backscattering  from  spherical  shell. 


BACKSCATTERING  FROM  SPHERICAL  SHELL 


Figure  2.  Contribution  of  mode  4  to  backscattering  from  spherical  shell. 
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The  second  example  involves  the  backscattering  off  the  end  of  a  capped  cylindrical 
shell.  The  ratio  of  the  mean  radius  a  to  the  length  L  of  the  shell  is  0.125,  and  the 
ratio  of  the  thickness  h  to  the  mean  radius  a  of  the  shell  is  0.01016  .  The  ends  of  the 
cylinder  are  capped  with  circular  disks  having  the  same  thickness  as  the  shell.  The 
material  properties  of  both  the  shell  and  endcaps  are 

•  Density  =  8977  Kg/m^ 

•  Young’s  modulus  =  2.09  x  10'*  N/m^ 

•  Poisson’s  ratio  =  0.308. 


The  material  properties  of  the  water  are 


•  Density  =  998  Kg/m^ 

•  Sound  Speed  =  1486  m/s. 


The  dotted  curves  in  figures  3  to  8  represent  the  contribution  of  selected  individual 
terms  in  the  summation  shown  in  equation  (19)  to  the  backscattering  off  the  end  of 
the  cylinder.  The  scattered  pressure  is  normalized  by  the  high-frequency  plane-wave 
approximation  to  the  rigid  backscattered  pressure.  The  plane  wave  approximation 
involves  setting  p,  =  pcv,  on  the  surface,  where  p,  and  v,  are  the  scattered  pressure 
and  normal  velocity  respectively.  For  a  rigid  body,  the  scattered  and  incident  normal 
velocities  on  the  surface  are  related  by  v,  =  — u,nc-  Thus,  the  plane  wave  approxima¬ 
tion  reduces  to  setting  p,  =  —pcvinc  on  the  surface.  The  far-field  scattered  pressure 
can  be  obtained  from  the  approximate  values  of  p^  and  v,  on  the  surface.  The  nor¬ 
malized  scattered  pressure  used  in  *his  report  is  very  similar  to  the  form  function, 
which  is  defined  to  be  the  ratio  of  the  scattered  pressure  to  the  geometric  acoustics 
approximation  of  the  scattered  pressure.  In  fact,  for  a  sphere,  the  two  normalizations 
are  the  same.  The  normalized  frequency  ka  is  defined  by  ka  =  2-Kfalc  where  /  is 
the  frequency,  a  is  the  mean  radius,  and  c  is  the  sound  speed  in  water.  The  dotted 
cu'^ves  in  figures  9  ctnu  10  represent  the  contribution  of  a  partial  sum  of  selected  terms 
in  equation  (19)  to  the  scattered  field.  The  solid  curve  in  each  of  these  figures  rep¬ 
resents  the  normalized  pressure  (in  dB)  with  the  rigid  scattering  excluded.  Figures 
11  to  16  show  the  in  vacuo  mode  shapes  of  the  modes  used  in  figures  3  to  8.  It  is 
clear  from  these  figures  that  certain  of  the  peaks  in  '  ’'.e  scattered  field  are  due  to  a 
single  in  vacuo  mode.  For  example,  „he  first  peak  at  about  ka  =  0.045  is  primarily 
due  to  mode  2  (see  figure  3).  As  can  be  seen  from  figure  11,  mode  2  is  predominately 
an  end  cap  mode.  However,  certain  features  such  as  the  hump  in  the  ka  =  d.3-0.5 
range  and  the  humps  in  the  ka  —  0.6-  1.0  range  require  the  contribution  of  several  in 
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vacuo  modes  (see  figures  9  and  10).  Figures  17  and  18  show  the  surface  pressure  and 
normal  velocity  on  the  cylinder  at  ka  —  0.04228  .  Notice  that  the  normal  velocity 
distribution  is  very  similar  to  the  in  vacuo  mode  shape  2.  Figures  19  and  20  show  the 
surface  pressure  and  normal  velocity  on  the  cylinder  at  ka  =  0.95136  .  Since  multiple 
modes  are  involved  at  this  frequency,  there  is  no  correspondence  between  the  normal 
velocity  distribution  and  any  single  mode  shape. 

Since  the  total  number  of  in  vacuo  modes  used  is  often  very  large,  it  is  useful  to 
first  scan  a  frequency  range  to  determine  which  modal  contributions  exceed  a  preset 
threshold.  The  results  of  a  sample  scan  are  shown  in  table  1  for  the  cylinder  problem, 
with  a  threshold  of  10  dfi  over  the  frequency  range  ka  =  0.025-0.3  . 
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Figure  4.  Contribution  of  mode  3  to  scattering  from  cylinder. 


9 


Figure  5.  Contribution  of  mode  5  to  scattering  from  cylinder 


Figure  6.  Contribution  of  mode  6  to  scattering  from  cylinder. 
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Figure  7.  Contribution  of  mode  14  to  scattering  from  cylinder. 


Figure  8.  Contribution  of  mode  15  to  scattering  from  cylinder. 
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Figure  9.  Contribution  of  modes  2,3,4,5,6,38,44. 


BACKSCATTERING  FROM  END  OF  CYLINDER 

ALL-RS .  PARm... 


Figure  10.  Contribution  of  modes  2,3,4,5,6,14,15,38,44. 
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MODES  5.  kas0.37$66 


Figure  14.  Cylinder  mode  shape  6. 
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MODE  =  14.  k»  =  3.0073 


Figure  15.  Cylinder  mode  shape  14. 


MODE=  15.  ki  =  3.1486 


Figure  16.  Cylinder  mode  shape  15. 
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CHIEF  VELDOTIES  FOR  kM  =  0i>4228 


Figure  17.  Normal  velocity  on  cylinder,  ka  =  0.04228. 


C(nEF  PRESSURES  FOR  k*  -  0.04228 
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Figure  18.  Surface  pressure  on  cylinder,  ka=0. 04228. 
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Figure  19.  Normal  velocity  on  cylinder,  ka  =  0.95136. 


CHIEF  PRESSURES  FOR  ka  =0.95136 
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CONCLUSIONS 


The  scattering  from  elastic  bodies  immersed  in  a  fluid  is  highly  variable  with 
frequency,  due  to  the  many  resonant  modes  of  the  structure.  Computational  models 
for  solving  structural-acoustic  scattering  problems  have  been  developed  at  NOSC.  The 
scattered  pressure  field  in  these  models  can  be  represented  as  a  sum  of  contributions 
from  in  vacuo  modes  of  the  structure.  Although  the  overall  spectral  response  of  a 
scatterer  is  typically  due  to  the  contribution  of  hundreds  of  such  modes,  it  has  been 
shown  that  the  response  in  narrow  frequency  bands  (especially  near  sharp  peaks)  is 
often  dominated  by  the  contributions  of  a  small  number  of  modes.  Auxiliary  programs 
have  been  developed  to  identify  and  display  these  modal  contributions  along  with  the 
associated  surface  pressure  and  velocity  distributions.  The  results  obtained  have 
provided  valuable  insight  into  the  mechanisms  responsible  for  features  observed  in 
the  farfield  response. 
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